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ABSTRACT 

We investigate the impact of the star formation efficiency in cluster-forming-cores (i.e. local SFE) 
on the evolution of the mass in star clusters over the age range 1-100 Myr, when star clusters un- 
dergo their infant weight-loss /mortality phase. Our model builds on the A-body simulation grid of 
iBaumgardt fc Kroupa I (|2007l ). Assuming a constant formation rate of gas-embedded clusters and a 
weak tidal field, we show that the ratio between the total mass in stars bound to the clusters over that 
age range and the total mass in stars initially formed in gas-embedded clusters is a strongly increasing 
function of the averaged local SFE, with little influence from any assumed core mass-radius relation. 
Our results suggest that, for young starbursts with estimated tidal field strength and known recent 
star formation history observed cluster-to-star mass ratios, once corrected for the undetected clusters, 
constitute promising probes of the local SFE, without the need of resorting to gas mass estimates. 

Similarly, the mass ratio of stars which remain in bound clusters at the end of the infant 
mortality/weight-loss phase (i.e. age > 50 Myr) depends sensitively on the mean local SFE, al- 
though the impacts of the width of the SFE distribution function and of the cor e mass-ra d ius re lation 
require more careful assessment in this case. Following the recent finding by iBastianl (|2008h that 
galaxies form, on the average, 8 % of their stars in bound clusters regardless of their star formation 
rate, we raise the hypothesis that star formation in the present-day Universe is characterized by a 
near-universal distribution for the local SFE. A related potential application of our model consists 
in tracing the evolution of the local SFE over cosmological lookback times by comparing the age 
distribution of the total mass in star clusters to that in field stars in galaxies where field stars can be 
resolved and age-dated. We describe model aspects which are still to be worked out before achieving 
this goal. 

Subject headings: galaxies: star clusters — galaxies: evolution 
tion — stellar dynamics 



galaxies: starburst — stars: forma- 



1. INTRODUCTION 

Star formation efficiencies (SFE) are not easy to come 
by observationally since they require reliable determina- 
tions of both the gaseous and stellar masses of a star 
forming region. Yet, the SFE is a fundamental parame- 
ter of star and star cluster formation processes and of the 
early evolution of star clusters. Local SFEs, i.e. the mass 
fraction of a giant molecular cloud core turned into stars, 
have been estimated for a handful of Gala ctic nearby gas- 
embedded clusters (jLada fc Lada I [2003) . They range 
from 10 to 30 per cents, which is an order of magnitude 
higher than what is inferred when the SFE is averaged 
over an entire giant molecular cloud (the so-called global 
SFE). Attempts to compare SFE estimates in the Milky 
Way and its neighbours, where they are best studied, to 
those in major merger-induced starbursts, involves the 
additional difficulty of scale differences. While for the 
former, SFEs are defined down to individual giant molec- 
ular clouds (GMC) and their dense star forming cores, 
only large-scale - often galaxy- wide - SFEs are measured 
for the latter owing to their large distances. Estimates of 
SFEs over up to kpc size nuclear regions in Ultra Lumi- 
nous IR Galaxies (=ULIRGs), all of which are late stages 

1 Belgian Science Policy Research Fellow and Humboldt Fellow 

2 Institute of Astrophysics & Geophysics, University of Liege, 
Allee du 6 Aout 17, B-4000 Liege, Belgium 

3 Argelander-Institut fur Astronomic, University of Bonn, Auf 
dem Hiigel 71, D-53121 Bonn, Germany 

4 Centre for Astrophysics Research, University of Hertfordshire, 
Hatfield ALIO 9AB, United Kingdom 



of massive gas-rich mergers, are reported to be 1 — 2 or- 
ders of magnit ude higher than global SFEs measured in 
spiral galaxies (|Gao fe Solomon 112004 their fig. 2b). 

Because the formation of compact and massive star 
clusters is observed to be closely associated to active 
star forming environments, prominent star formation 
episodes in the history of galaxies are best traced by the 
age distribution of the i r star c lusters (IWest et ai~l 20041 : 
Fritze-v. Alvensleben I 120041: iBrodie fc Strader I 120061 : 



Kotulla. Fritze fc Anders 1120081) . Although the vast ma- 



jority of stars actually form in gas-embedded star clus- 
ters, galaxies are predominantly made of field stars. In 
this respect, a peak in a star cluster age distribution only 
tells us when an epoch of enhanced star cluster forma- 
tion took place, not yet how much total stellar mass was 
involved. To derive the total mass in stars formed in a 
starburst older than ~ 3—4 Gyr from the integrated spec- 
tru m of a galaxy is, in the present state-o f-the-art, hope- 
less (jFritze-v. Alvensleben fc Lilly I [2007). Knowing the 
mass fraction of all star formation still in bound clusters 
as a function of star cluster age, by properly accounting 
for cluster destruction and mass loss, would allow us to 
carry the analysis of galaxies'evolution one step further 
and to estimate the history of their overall star formation 
rate. We now have a fairly mature view of how to deci- 
pher the rates of bound cluster formation and evolution of 
a given star cluster system, in es sence based on the clus- 
ter a ge and mass distributions ijParmentier fc de Griis I 
2008). That is, we are able to assess the amount of field 
stars arising from the secular evolution (i.e. the prefer- 
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ential depletion of low-mass stars due to tidal stripping 
and internal two-body relaxation) of star clusters which 
survived the expuls ion of their re s idual star forming gas. 
As pointed out by iMcLaughlin I (|1999( ), however, secu- 
lar evolution does not necessarily affect the total mass in 
star clusters significantly. As for the star clusters of the 
Galactic Old Halo and of the Large Magellanic Cloud for 
instance, the observed total mass in clusters is less than 
an order of magnitude smaller than their total mass at 
th e onset of secular evolution (see table 3 and table 2 
in iParmentier fc Gilmorel 120051 : iParmentier fc de Griis I 
120081 respectively). If the vast majority of stars form in 
gas-embedded star clusters, this implies that the bulk of 
field stars is given off at an earlier stage, very likely as 
a result of cluster dynamical response to the expulsion 
of the residual star forming gas out of their parent cores 
due to massive star activity. Because of its resulting shal- 
lower potential, an exposed cluster attempts to reach a 
new equilibrium, which results into the loss of stars ( "in- 
fant weight-loss") and, when the local SFE is too low, 
the eventual disru ption of the cluster ("infant mortal- 
ity") (see Fig. 1 in iGoodwin fc Bastian Il2006t I Goodwin! 
2008). This violent relaxation phase, which lasts for at 
most 50Myr, appears therefore as the prime driver of 
the ratio of the total mass in bound clusters to the total 
mass in stars. In what follows, we assume that all stars 
form in gas-embedded clusters, i.e. we assume that there 
is no mode of field star formation or loose stellar asso- 
ciation formation. We will further comment about this 
hypthesis in our discussion (section 3). 

Young starbursts are the best places to infer the im- 
pact of cluster violent relaxation on the cluster-to-star 
mass ratio. Based on HST UV imaging of a sample 
of starburst galaxies, the morphology of which ranges 
from blue compac t dwarfs to ultra-lum inous merging far- 
infrared galaxies. iMeurer et al. I (|1995f ) find that, on av- 
erage, about 20 p er cents of the UV lu minosity comes 
from star clusters. Ide Griis et al. I (|2003l ) study the Mice 
and Tadpole interacting systems by means of a pixel-by- 
pixel analysis of their color-magnitude and color-color di- 
agrams. In the Mice interacting galaxies, they find that 
star clusters are predominantly associated with regions 
of active star formation, where they contribute > 40 per 
cents of the total flux, irrespective of wavelength (but see 
our discussion section for caveats regarding these values). 
Similar flux ratios are inferred for the Tadpole system, an 
example of a galaxy encounter between two unequal-mass 
galaxies. These high luminosity fractions, presumably 
representative of large mass fractions, demonstrate that 
star clusters contribute significantly to the total stellar 
light given off by young starbursts. While being higher 
than the initial Galactic halo globular cluster mass frac- 
tion as compared to halo field stars, i.e. ~ 0.1, these ra- 
tios are also conspicuously smaller than unity, thus con- 
firming that the dissolution of star clusters takes place 
on a short time-scale. 

Analyses of star cluster systems on the basis of multi- 
band imaging are feasible out to Virgo cluster distances 
and have the potential to provide a comprehensive view 
of the individual star formation, chemical enrichment 
and mass assembly histories of their parent galaxies over 
the past Hubble-Time, provided that the last missing 
issue, establishing the ratio of star formation ending 
up in bound star clusters, is solved. This is what we 



attempt here, that is, we investigate what is/are the 
driving parameter(s) behind the mass fraction of clus- 
ters as compared to field stars. Specifically, we build 
on the A^-body simulation g rids recently obtained by 
Baumgar dt fc Kroupa I (|2007f ) to compute the evolution 
over the age range 1-100 Myr of the mass in clusters and 
to compare it with the overall star formation rate arising 
from the same star formation episode. We make predic- 
tions about the cluster formation efficiency (i.e. the ratio 
between the total mass in clusters over the age range of 
relevance and the total mass in star forming gas cores) 
and the ratio between the total mass in clusters and the 
total mass in stars. And we highlight their relation to 
the local SFE. 

2. METHOD AND RESULTS 

The starting point of our simulations is an ensemble of 
molecular cores, the high-density regions of GMCs asso- 
ciated to star formation. The core population is statisti- 
cally described in the following way. The core mass spec- 
trum follows a power-law dN oc m a dm of spectral index 
a = — 2 and of lower and upper mass limits m; = 100 Mq 
and m up = 10 7 M Q , respectively. Each cluster forming 
out of one of these cores is assigned an age t randomly 
drawn from a uniform age distribution ranging between 1 
and 100 Myr. Cluster ages will be used to compute their 
mass and to build the temporal evolution of the mass in 
star clusters through the violent relaxation phase. Each 
of our Monte-Carlo simulations encompasses 10 6 cores 
with mass m c , radius r c (the core radius is defined be- 
low) and corresponding cluster age t. Combined to the 
assumed core mass spectrum, this core number is equiv- 
alent to a total mass of star forming gas of 1.1 x 10 9 Mq. 
We note that the core number has been chosen so as 
to limit Poissonian noise in the gas-embedded cluster 
formation rate we will derive, rather than as a match 
to a particular galaxy. We account for star formation 
with a time-independent Gaussian probability distribu- 
tion function of the SFE of mean e and standard devia- 
tion <7 e , which we denote P(e) = G(e,a e ). We emphasize 
that e is the local SFE, that is, the core mass fraction 
turned into stars, as opposed to global SFE measured 
on galactic-wide or GMC-wide scales. The correspond- 
ing constant star formation rate is thus ex 11 MQ.yr^ 1 , 
equivalent here to the total stellar mass formed per unit 
of time in gas-embedded clusters. 

To obtain the stellar mass m c \(t) still bound to a clus- 
ter of age t through the infant mortality/infant weight- 
loss phase, the gas-embedded cluster mass e x m c is mul- 
tiplied by the bound fraction Fb of stars at age t: 

m c i(t) = F h (t) e m c . (1) 

Stellar groups of which the stellar mass is less than 
100 Mq are not considered as star clusters and their mass 
is not accounted for in the total cluster mass. Because 
of their small total mass, this is of negligible impact 
on the clust er-to-star mass ratios we deriv e. In the N- 
body grid of lBaumgardt fc Kroupa I (|2007l ). from which 
our star cluster system simulations are derived, Fb(t) is 
defined as the stellar mass fraction of the initially gas- 
embedded cluster which resides within the instantaneous 
cluster tidal radius. It depends on the local SFE e, on 
the age t of the cluster, on the gas removal time-scale 
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tgr/tctoss (expressed in unit of a protocluster crossing- 
time) and on the strength of the external tidal field. The 
older the cluster and/or the lower the star formation ef- 
ficiency and/or the quicker the gas expulsion and/or the 
stronger the host galaxy tidal field, the lower the bound 
fraction Ft,. To recover this four- fold dependence of Ft,, 
we interpolate the output files of the TV-body sim ulation 
runs carried out bv lBaumgardt fc Kroupa I (|2007l ). Their 
grid of results provides the evolution with time of Ft, fol- 
lowing gas expulsion, as a function of e, of TGR,/T cross 
and of the ratio of the half-mass radius to the tidal ra- 
dius of the gas-embedded cluster rh/rt , which quantifies 
the strength of the tidal field. 

The local SFE and the age t are sampled from 
their respective distribution functions. The gas re- 
moval time-scale tor. /7v- rr ,sK is defined by equation 6 in 
iParmentier et al~l (|2008f ): it depends on the local SFE 
e and on the core mass m c and radius r c . It remains 
unclear how r c scales with m c . Theoretical expectations 
are that the mass-radius relation of virialized cores obeys 
r c oc ■m L J' 1 (e.g. lHarris fc Pudritz I [l~994f ) . However, this 
is at variance with observations that the radii of gas- 
embedded clusters and recently exposed clusters do not 
vary significantly (r ~ 1 pc) over severa l orders of mag - 
nitude of cluster masses (see Table 1 in iKroupa 1120051) . 
Considering how debated this issue remains, in what fol- 
lows, we consider three distinct cases: cores with con- 
stant radius (r c oc mP c ), constant density (r c oc m^ 3 ) 

l li 

and constant surface density (r c oc m c , i.e. virialized 
cores) . Corresponding simulation results are respectively 
shown as cases (a), (b) and (c) in Figs [TJ [5] and [3J and in 
Table Q] As for the normalizati on of these relations, we 
refer to Figure I of iTan I (|2007f ). which shows the mass, 
radius, volumic mass density and surface mass density 
of gas-embedded Galactic protoclusters. In case (a), we 
assume that the protocluster radius is a good proxy to 
the core radius and we adopt r c — 0.7pc. In case (b), 
we fit by eye a constant volumic mass density line to the 
Galactic protoclusters and we multiply the correspond- 
ing density by about a factor of 3, so as to account for 
the mass contribution of the residual star forming gas 
(i.e. we assume a mean local SFE of 0.33). This gives: 
r c = 0.026 x (mc/IAl©) 1 / 3 pc. In case (c), adopting a 
similar approach for a constant surface mass density, we 
find r c = 0.008 x {raJlMQ) 1 ' 2 pc. Lower normalizations 
(i.e. smaller radius for a given mass) would slow down 
gas expulsion, thereby increasing the fraction of stars re- 
mai ning bound to the cluster (see equation 6 and Figure 
1 in IParmentier et al. 1120081 ). In order to highlight the 
impact of gas expulsion only on the cluster age distribu- 
tion, we limit our simulations to the case of a weak tidal 
field. The vast majority of our gas-embedded clusters 
have Thjrt < 0.03. Those characterised by rh/rt > 0.03 
(which appear only when r c = 0.7 pc; case (a)) have a 
mass of order 100 Mq and become therefore rapidly ir- 
relevant to the cluster age distribution. As an example, 
a cluster with a mass of 10 4 Mq and a half-mass radius 
of 0.8 pc orbiting the Milky Way halo at a galactocentric 
distance of 37kpc has rh/rt ~ 0.01. In a forthcoming 
paper, we will extend our model to stronger tidal fields, 
so as to assess the impact of the environment on the early 
star cluster history. 

Once the bound fraction Ft,(t, e, tgr/t ci . oss , rh/rt) has 



been estimated, the instantaneous (bound) star cluster 
mass corresponding to each star forming core is obtained 
via equation [1] The corresponding temporal evolution 
of the mass in star clusters is then straightforwardly 
built. We have simulated 18 different star cluster sys- 
tems, combining the three core mass-radius relations de- 
scribed above with two mean SFEs (e = 0.25 and 0.40) 
and three standard deviations (cr e = 0.01,0.04,0.07) for 
the SFE distribution function P(e). Results of these sim- 
ulations are illustrated in Figs. [T] and [D 

Each panel shows the time evolution of the mass 
in dense core gas available to star formation per time 
unit (solid lines with asterisks), the star formation rate 
(SFR; dash-dotted lines with plain circles), and the re- 
sulting age distribution of the mass in star clusters 
('SCs': dotted lines with open triangles, squares, cir- 
cles for <7 e = 0.01,0.04,0.07, respectively). All mass 
counts are in unit of 1 MQ.yr . Because of model hy- 
potheses, histories of the SFR and of the core forma- 
tion rate are time-independent. In contrast, the com- 
bined infant mortality and infant weight-loss of star clus- 
ters is clearly highlighted. About 8Myr after gas ex- 
pulsion, the star cluster age distribution starts decreas- 
ing as a result of unbound stars reaching cluster tidal 
limits. At an age of ~ 30-50 Myr, violent relaxation is 
over and the age distribution stabilizes (provided that a 
fraction of the initial gas-embedded cluster population 
survives). That the age distribution eventually reaches 
an almost constant value is also due to our assump- 
tion of a weak tidal field for the simulated star clus- 
ter systems, which is appropri ate for, e.g, the Small 
and Large Magellanic Clouds (seelBoutloukos fe Lamers I 
120031 IParmentier fc de Griis 1 120081 . respectively). In 
dense environments, such as the Whirlpool galaxy M51 
for which the dissolution time-scale - due to secular evo - 
lution - of a 10 4 M Q is ~ 100 Myr (jGieles et al. Il2005| ). 
the cluster age distribution will keep declining, albeit at 
a different rate, as a result of the sole cluster secular evo- 
lution (two-body relaxation and tidal interactions with 
the host galaxy). 

Figs. [T] and [2] show that the post-violent-relaxation 
cluster age distribution is a sensitive function of P(e) — 
G(e, cr e ). (Note that the lower limit of the y-scale in 
Fig. [2] is much higher than in Fig. Q]). In case of instan- 
taneous gas expulsion (tqr << r cross ) and weak tidal 
field, the SFE threshold eth required to get a bound 
cluster at an age of > 50 M yr is e t h — 0.33 (Fig. 1 in 
IParmentier fe Gilmore 1120071 ) . Therefore, any P(e) func- 
tion leading to a negligible contribution of the e > e t h 
regime (e.g. P(e) = (7(0.25,0.01)) leads to the disrup- 
tion of the whole cluster population on a time-scale < 
30 Myr. In contrast, when e = 0.40, the vast majority of 
the clusters survive, albeit depleted by about 65 per cents 
of their mass, i.e. F b (t = 50 My r, e = 0.40) ~ 0.35 (Fig. 1 
in lParmentier fc Gilmore ll2007h . The decrease with time 
in the mass in star clusters is, in this case, an imprint 
of infant weight-loss only. Generally speaking, when a is 
small, the bound fraction Ft,(t = 50 Myr, e) is an excel- 
lent proxy to the ratio between the age distributions of 
the masses in clusters and in stars for t > 50 Myr. When 
e < eth, the final cluster-to-star mass ratio depends on 
a since it determines how large the contribution of the 
e > eth regime is (see also the discussion in the next 
section) . 
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Fig. 1. — Time evolution of the mass in clusters over the course of 
violent relaxation. The mass of core gas available to cluster forma- 
tion per time unit (solid lines with asterisks) and the distribution 
function of the star formation efficiency are assumed constant in 
time, resulting into a constant star formation rate (dash-dotted 
lines with plain circles). The age distribution of clusters in terms 
of mass counts is illustrated for different core mass-radius rela- 
tions (panels a, b and c) and for different standart deviations of 
the star formation efficiency probability distribution (dotted lines 
with open circles, squares and triangles). The mean star formation 
efficiency is e = 0.25 



The assumption of either core mass-radius relation 
mostly impacts on the history of the mass in clusters 
when the SFE is low. Considering the case of P(e) = 
G(0.25, 0.04) (dotted curves with open squares in FigQ]), 
the age distribution at t > 50 Myr is increasing when 
the core mass-radius relation gets shallower (from panel 
c to a; in the extreme case, when the core radius is con- 
stant). The origin of this behaviour is to be sought in 
Fig. [3] which depicts the relation between the gas expul- 
sion time-scale tgr/t ctoss and the core mass m c for the 
different core mass-radius relations used in our simula- 
tions. The thick horizontal line at tgr/t c10SS ~ 1/3 parts 
the regime of explosive gas expulsion from the adiabatic 
one. A ctually, inspection of Fig. 1 in iParmentier et al71 
(2008), which shows Ff,(t — 50 Myr, e, TGR,/T cr0 ss) for a 
weak tidal field, illustrates that as long as TGR/T croS s < 
1/3, the bound fraction is constant and characteristic of 
instantaneous gas removal. Once the protocluster gas is 
expelled on a time-scale slower than tgr/t c10SS ~ 1/3, 
the bound fraction is steadily increasing with the gas ex- 
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pulsion time-scale. Fig. [3] shows that when : 
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Fig. 2. — Same as in Fig.[T] for a mean star formation efficiency 
of e = 0.40. Note the difference in the y-scaling with Fig. [T] 



gas cores undergo explosive gas expulsion. In contrast, 
when r c is constant, cores more massive than 10 5 M Q 
undergo adiabatic gas expulsion 5 . The bound fraction, 
for any given age, is then higher than for instantaneous 
gas expulsion and is a sensitive function of tgr/tctosb- 
As a result, the stellar mass fraction remaining in bound 
clusters gets higher as the contribution of the adiabatic 
regime is growing. 

A puzzling feature in the history of the cluster mass 
when e = 0.25 is that it increases over the age range 50- 
100 Myr for shallow core mass-radius relations (e.g. open 
triangles in panel a and open squares in panels a and b of 
Fig- [J)! while one would expect a stabilizing age distri- 
bution. This is a consequence of the size-of-sample effect 
for the core mass spectrum, combined to the ability of 
high- mass cores to undergo adiabatic gas expulsion when 

r c oc or when r c is constant (Fig. [3]). Let us consider 
the oldest age bin for a — 0.01 and r c = constant (right- 
most open triangle in panel a). It covers the age range 
7.8 < log(t) < 8.0 and is contributed to by three clusters 
only, all arising from cores more massive than 5 x 10 6 M . 
Because of the assumed constant core radius, these very 

5 It is worth keeping in mind that the quantitative results of 
Fig. [3] depend on the normalisations we have adopted for each of 
the core mass-radius relations. A core radius twice as large for 
any given core mass would make gas ex pulsion almost twice as fast 
(equation 6 in Parmentier ct al. (2008)) 
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Fig. 3. — Relations between the gas expulsion time-scale in unit 
of one protocluster crossing time and the cluster forming core mass, 
for different core mass-radius relations. The steeper the core-mass 
radius relation, the weaker the dependence of the gas expulsion 
time-scale on the core mass. A sample only of the one million 
cores simulated is shown 



massive cores undergo slow gas expulsion and, therefore, 
retain a non-zero fraction of their stars (of order ~ 50 per 
cents for two of them) despite an SFE e ~ 0.25. Less mas- 
sive cores, owing to their shallower gravitational poten- 
tial, result into disrupted clusters (Fb = 0). At younger 
age (7.4 < log(t) < 7.6), constant logarithmic age bins 
cover shorter linear age ranges, resulting into a smaller 
number of cores per bin, and a decreased probability of 
including very massive cores able to provide bound star 
clusters. 

3. DISCUSSION 

The previous section has highlighted the steep rela- 
tionship binding the mean local SFE e and the history 
of the mass in star clusters under the assumptions of a 
weak external tidal field and of a constant gas-embedded 
cluster formation rate. Table Q] lists, for the 18 runs of 
Monte-Carlo simulations, the global (i.e. integrated over 
the age range 1-100 Myr) cluster formation efficiency e c i 
and cluster-to-star mass ratio M c i/M st . The cluster for- 
mation efficiency e c / is defined as the ratio between the 
total mass in stars bound to the clusters over the age 
range 1-100 Myr and the total mass in core gas. These ra- 
tios appear robust with respect to both the assumed core 
mass-radius relation and the standard deviation <r e of the 
SFE distribution function. Tabled demonstrates that, for 
this set of simulations, the cluster formation efficiency 
e c i and the cluster mass-to-star mass ratio M c i/M st are 
sensitive tracers of the local mean SFE e. Depending on 
the core mass-radius relation and SFE standard devia- 
tion, the ratio between the total mass in clusters and 
that in stars is 4 to 7 times higher when the mean SFE 
is e = 0.40 than when e = 0.25. This tight relation is a 
direct consequence of residual star forming gas expulsion 
out of protoclusters. For a population of clusters which 
have survived their violent relaxation and which are not 
yet affected by secular evolution significantly (or equiva- 
lently which are subject to a weak tidal field), the ratio 
between the total mass in clusters and the total mass in 



TABLE 1 

Ratio between the total mass in clusters and the total 
mass in stars (m c i/m a t) and cluster formation efficiency 
€ c i (i.e. ratio between the total mass in clusters and the 

total mass in core gas j for all simulations shown in 
Figs. [T]and [2] All total masses are integrated over the 
age range 1-100 Myr . A weak tidal field and a constant 

GAS-EMBEDDED CLUSTER FORMATION RATE ARE ASSUMED. 
BRACKETTED VALUES CORRESPOND TO CLUSTER-TO-STAR MASS 
RATIOS FOR CLUSTERS BRIGHTER THAN My = —7.3, THEREBY 
EMULATING A DETECTION LIMIT 







e = 0.25 
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0.01 


0.04 


0.07 


0.01 


0.04 


0.07 






(a) r c 


= constant 








M cl /M st 


0.06 


0.06 


0.10 


0.42 


0.42 


0.43 




(0.04) 


(0.04) 


(0.07) 


(0.27) 


(0.27) 


(0.27) 


e c i 


0.02 


0.02 


0.03 


0.17 


0.17 


0.17 
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1/3 

r c oc m c 








M cl /M 3t 


0.07 


0.08 


0.10 


0.36 


0.36 


0.39 




(0.06) 


(0.06) 


(0.07) 


(0.23) 


(0.22) 


(0.25) 


ed 


0.02 


0.02 


0.03 


0.15 


0.14 


0.16 






(c) 


1/2 
r c oc m c 








Mci/Mst 


0.08 


0.08 


0.10 


0.36 


0.36 


0.40 




(0.06) 


(0.07) 


(0.08) 


(0.23) 


(0.24) 


(0.27) 


£d 


0.02 


0.02 


0.03 


0.15 


0.15 


0.16 



stars is simply given by: 
— (50 < t < 100 Myr) 



Fbound X f X M e 

e x M c 



Abound : 



(2) 

where M c is the total mass in core gas, Fb oun d and e are 
the bound fraction F 6 (50 < t < 100 Myr, e) and the SFE 
averaged over the whole population of clusters, respec- 
tively. The cluster mass fraction, as compared to the 
total mass in stars formed in the same star formation 
episode, is thus governed by the post-violent-relaxation 
bound fraction which is itself a steep function of the 
local SFE (see Fig. 1 of IParmentier fc Gilmore 1 120071 
for the case of a weak tidal field) 6 . Measurements of 
the flux ratio from star clusters relative to field stars 
in young starbursts combined to improved models ac- 
counting, when required, for a significant external tidal 
field and a non-constant gas-embedded cluster formation 
rate thus constitute a promising way of probing the lo- 
cal SFE in active star forming environments, without 
resorting to gas mass estimates. Prerequisites include 
that the age range of the star cluster system includes the 
post-violent-relaxation age range, i.e. [50,100] Myr for a 
weak tidal field. If the tidal field is stronger, cluster vio- 
lent rela xation is faster since protoclu ster tidal radii are 
smaller ( Baumgardt &; Kroupa Il2007h . 

Table [T] also shows that cluster mass fractions as high 
as 40 per cents are not surprising: they are the signature 
of a local SFE only slightly in excess of the threshold 
£th required to form a bound star cluster. If e ~ 0.50, 
the cluster mass fraction leaps to 70 per cents. It is 
even higher in case of slow gas expulsion (Fig. 1 in 
IParmentier fc Gilmorel l2007f ) . as is the case for small 

6 Note that M c i/M st differs from in TableQ]even for e = 0.25, 
cr £ = 0.01 and explosive gas expulsion because the cluster-to-star 
mass ratio in this table is integrated over the age range 1-100 Myr 
and, therefore, accounts for clusters still little affected by violent 
relaxation 
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Fig. 4. — Age-mass diagram of the fiducial star cluster system 
with r c oc ml /3 and P(e) = G(0.40, 0.04). The solid line corre- 
sponds to a constant magnitude of My 111 = —7.3. A fraction only 
of the clusters arising from the one million cores simulated is shown 



size massive cores. 

However, it must be kept in mind that these cluster 
mass fractions are upper limits. Firstly, at a given 
age, the presence of a strong tidal field, as may be 
expected in dense star formin g environments, would 
make them lower (see Table 1 in iBaumgardt fc Kroupal 
|2007t ). Secondly, the census of star clusters in a distant 
galaxy is necessarily incomplete and the detection limit 
implies that the lower cluster mass range is missed. 
We now investigate this second point, that is, how 
much the cluster-to-star mass ratio is affected by the 
presence of a detection limit. Should it not be corrected 
for, the M c i/M st ratio will be underestimated and an 
erroneously low local SFE estimate will be inferred. We 
illustrate this issue with a 100 per cent-completeness 
limit at M v m = —7.3, which is compa rable to what 
has been found by iBastian et al. I ([2005) in their HST 
study of the M51 star clusters, at a distance of 8.4 Mpc. 
All star cluster masses were converted into luminosities 
with the GALEV model (jSchulz et al. Il2002f ) and an ad- 
ditional cluster-to-star mass ratio was computed, which 
takes into account clusters brighter than My m = —7.3 
only. These are the bracketted values in Table [TJ As 
an example, the age-mass diagram of the star cluster 

1 /3 

system with r c oc m c and P(e) = G(0.40, 0.04), along 
with the adopted detection limit, is shown in Fig. 0J 
Although the vast majority of the clusters are fainter 
than M v m = -7.3, the M c i/M st ratio is reduced by 
about one third only. This is due to the choice of 
a power-law of spectral index —2 for the core mass 
spectrum. While clusters with a mass in the range 
of, say, 100 Af© < m c i < 3000 Af© (where the upper 
limit is chosen as a rough averaged estimate of the 
cluster mass at the detection limit) represent 97 per 
cents of the cluster number, they contribute by 30 per 
cents only to the total mass of the star cluster system. 
Their non-detection therefore does not severely impact 
on the M c i/M s t ratio. The effect will of course be 
more prominent as the detection limit gets brighter. 
Knowledge of the detection limit and of the age range 
covered by the star cluster system, combined to detailed 



simulations such as those carried out in this paper, are 
necessary to convert the observed M c i/M st ratio into 
the actual one. We also note that in distant galaxies, 
observed cluster-to-star mass ratios may be affected by 
confusion effects, i.e. groups of otherwise unbound stars 
being misinterpreted as genuine bound star clusters, 
which results in an overestimated total mass in clusters. 
How these two effects - detection limit and confusion 
- affect the c l uster- to-st ar mass/flux ra t ios de rived by 
iMeurer et al~1 (|1995[ ) and lde Griis et al. I (|2003f ) remains 
to be seen. 

In a recent study, IBastian 1 (|2008f) argues that the for- 
mation of bound clusters - those which survive violent 
relaxation - represents, on the average, 8 % of the to- 
tal star formation of a galaxy, irrespective of its SFR. 
// the star clusters displayed in the relation between 
the present SFR of gala xies and the ma gnitude of their 
brightest young clusters (|Bastian 1120081 his Fig. 1) have 
actually reached the endpoint of their violent relaxation, 
then we can use this result to estimate the correspond- 
ing local SFE. If all stars form in gas-embedded clus- 
ters, the 8 % ratio then corresponds to our bound frac- 
tion Ffjound, which gives a back-of-the-envelope estimate 
for the local star formation efficiency of e ~ 35 % (assum- 
in g explosive gas expulsion and weak tidal field; Fig. 1 
in iParmentier fc Gilmorel 120071), rega r dless of the host 
galaxy. According to lLada fc Ladal ((2003), however, 
gas-embedded clusters account for 70-90% of all stars 
formed in Galactic disk GMCs. If we now assume that 
30 % of all star formation in a galaxy takes place in the 
field or within loose stellar associations, the bound frac- 
tion, i.e. the ratio between the total mass in bound star 
clusters and that in gas-embbeded star clusters, becomes 
Fbound — 8%/0.7 ~ 11% which, by virtue of the sharp 
dependence of Fi, oun d on e (especially in the e > e t h 
regime), leads to an only marginally higher mean local 
SFE. 

The combination of the result of iBastianl (|2008| ) with 
our model provides therefore a hint for a local SFE 
which does not vary significantly among galaxies of 

the local Univer se. This is consistent with fig. 2a in 

iGao fc Solomon 1 (|2004fl which shows that the ratio be- 
tween the total far-IR luminosity Lib,-, a tracer of the 
SFR, and the HCN line luminosity Lhcn, a tracer of 
the dense star forming gas, is roughly constant (although 
with much scatter) over 4 orders of magnitude in SFR, 
from almost passive spirals to the strongest ULIRGs. It 
is worth noting that the local SFE estimate quoted above 
is only slightly higher than the SFE threshold tth for 
bound cluster formation in the regime of explosive gas 
expulsion. Therefore, even slight variations around the 
average value e ~ 0.35 are conducive to significant vari- 
ations in the bound cluster-to-star mass ratio (say, from 
% to 35 % for an SFE range of e <30 % to e -40 %). 

Merging galaxies have high SFRs, the origin of which 
is in the large fraction of molecular gas "trapped" in 
the dense phase of star forming cores, rather than in 
the more diffuse phase of GMCs. For instance, ULIRGs 
have conspicuously high mass ratio of dense star form- 
ing gas compared to their total molecular g traced 
by the HCN line luminosity to CO luminosity ratio 
Lhcn /Leo (|Solomon et all Il997t ISanders fc Mirabel I 
1996). This is the very reason of their l-to-2 orders 
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of magnitude higher SFR compared to spiral galax- 
ies (|Gao fe Solomon 1 12004 their fig. 4). Because of the 
observed correlation between the SFR and the vi sual ab- 
solute magnitude of the brightest young cluster (jLarsen I 
l2002tlWeidner. Kroupa fc Larsen l)200l . increased SFRs 
in turn gives rise to the formation of more massive clus- 
ters, which are more likely to survive several Gyr of secu- 
lar evolution. This is why major star formation episodes 
of a galaxy are well-traced by peaks in its cluster age 
distribution. 

4. MODEL CAVEATS 

Table Q] shows, for the models explored in section 2 
(i.e. weak tidal field and constant gas-embedded cluster 
formation rate), the ratio between the total mass in clus- 
ters and the total mass in stars formed in gas-embedded 
clusters over the age range 1 to 100 Myr (i.e. this ratio 
accounts for clusters of which the stellar content is still 
unaffected by violent relaxation), ft appears that these 
cluster-to-star mass ratios do not sensitively depend on 
either the standard deviation er e of the SFE distribution 
function or the core mass-radius relation. In contrast, the 
ratio between the total mass in bound clusters and that 
in stars (i.e. the cluster-to-star mass ratio integrated over 
the age range 50-100 Myr) does depend on the fine details 
and underlying model assumptions, as can be seen from 
Fig. Q] and Fig. [2] Because an increased bound cluster- 
to-star mass ratio can be caused by either an increased 
mean local SFE e, an increased SFE distribution function 
standard deviation cr € , or a shallower core mass-radius re- 
lation m c — r c , the model is prone to degeneracies. As for 
the mass-radius relation, the influence of gas expulsion 
on the slope of the cluster mass-radius relation as com- 
pared to that of their progenitor cores must be studied. 
In a forthcomin g paper, we will apply our model to the 
data of Fig. 1 of iBastianl (|2008f ) so as to constrain more 
tightly the distribution function of the local SFE. That 
is, we will investigate in greater detail what is the range 
in local SFE permitted by the data scatter around the 
linear relation between the logarithm of the present SFR 
and the visual absolute magnitude of the brightest young 
cluster in galaxies. 

We also remind the reader that our present re- 
sults are valid only in the case of a weak tidal field. 
Strong tidal fields both decrease the ability of gas- 
embedded clusters to survive infant weight-loss and 
shorten the violent relaxation phase, in essence be- 
cause a smal ler cluster tidal radius enha nces cluster 
loss of stars ()Baumgardt fc Kroupa I [20071 ). One way 
of assessing the tidal field strength in a given (re- 
gion of a) galaxy is to estimate the cluster dissolution 
time-scale due to secular evolution, based on the age 
and ma ss distributions of clusters older than 50 Myr 
(see e.g.lBoutloukos fc Lamers I l2l)0l iGieles et al. 1 120051: 
IParmentier fe de Grijs 1120081 ). One can then match this 
time-scale to the external tid al field via equation 7 of 
iBaumgardt fc Makino 1 (|2003l) . 

Finally, our simulations assume that the mass of dense 
core gas available to star formation per unit of time and 
the SFR remain constant with time. If a galaxy is en- 
tering a ULIRG phase, the amount of dense star forming 
gas will increase as time goes by, boosting the SFR and 
CFR at young ages. Should this effect not be corrected 
for, one could - wrongly - infer that the decrease in the 



cluster age distribution towards old ages is entirely due 
to cluster destruction and, therefore, underestimate the 
local SFE. 

5. CONCLUSIONS 

Building on the iV-body simula t ion g rid recently ob- 
tained by Bau mgardt fc Kroupa I (|2007[ L we have per- 
formed detailed simulations of the early (i.e. up to an age 
of 100 Myr) evolution of star cluster systems. In this first 
contribution, we have adopted the following assumptions: 
star forming cores with a power-law mass spectrum of 
slope —2, constant gas-embedded cluster formation rate, 
Gaussian distribution function of the local SFE and weak 
external tidal field. Besides, star formation is assumed to 
mostly take place in gas-embedded clusters. Under these 
assumptions, we have highlighted the strong impact of 
the local SFE on the temporal evolution of the mass in 
clusters and on the ratio between the total mass in clus- 
ters and the total mass in stars. Within less than 10 Myr 
after gas expulsion, the age distribution of the mass in 
clusters decreases as a result of cluster violent relaxation 
and the amplitude of the decline depends strongly on 
the adopted SFE distribution function. This is because 
the fraction of stars remaining bound to a cluster af- 
ter its violent relaxation is shar ply increasing with the 
SFE e once e > 0.33 (Fig. 1 in IParmentier fc Gilmorel 
l2007f ). The ratio M d /M st between the total mass m 
stars still bound to the clusters and the total mass in 
stars formed in gas-embedded clusters over the age range 
of 1 — 100 Myr constitutes a sensitive probe of the local 
SFE which characterized their progenitor cores. Besides, 
this ratio is only weakly affected by the choice of a core 
mass-radius relation. We caution however that (i) ob- 
servationally derived cluster-to-star mass ratios account 
for detected clusters only; (ii) ratios displayed in Table 
[I] will be different for other tidal field strengths and non- 
constant gas-embedded cluster formation history. Mea- 
suring the flux ratio given off by star clusters in a young 
starburst and converting it in a cluster-to-star mass ra- 
tio therefore constitutes a promising way of inferring the 
local SFE, provided that the observed M c i/M st is turned 
into the actual one and that time-varying gas-embedded 
cluster formation rate and stronger tidal fields are - if 
necessary - accounted for via generalized versions of the 
present model. 

Another potential application of our model is to trace 
the evolution of the local SFE over cosmological look- 
back times. This can be done, for galaxies close enough 
so that their field stars can be resolved and age-dated, 
by comparing the history of the bound cluster forma- 
tion rate (i.e. the age distribution of the cluster mass at 
the end of violent relaxation, obtained by correcting the 
observed age distribution for se cular evolutionary mass 
losses and fading effects, see e.g. IParmentier fc de Griis I 
(2008)) on the one hand and the history of the SFR on 
the other hand. This bound cluster-to-star mass ratio 
is our bound fraction Fb ounc [ (see equation [2}, equivalent 
to model outputs over the age range 50-100 Myr. This 
application, however, requires tighter constraints on the 
slope of the core mass-radius relation, the width of the 
local SFE distribution and an estimate of the external 
tidal field strength. 

Combined with the recent finding of IBastianl {2008) 
that, on the average, galaxies form 8 % of their stars in 
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bound star clusters, the relation between the bound frac- 
tion Fbound and the local SFE e constitutes a hint towards 
an almost universal local SFE in the present-day Uni- 
verse. In a forthcoming paper, we will apply our model 
to the observed relation between the present-day SFR 
of galaxies and the absolute visual magnitude of their 
brightest young clusters so as to shed better constraints 
on the present-day local SFE distribution function. 
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